On a new twist in the dynamics of Bose-Einstein condensation 
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The dynamical process of Bose-Einstein (BE) condensation of a dilute ultracold gas is found to 
occur in two distinctly different forms. If the particle flux into the low-energy modes is sufficiently 
strong, the Bose gas behaves predominantly like an incompressible fluid. It exhibits chaotic (Vinen) 
superfluid turbulence which marks the presence of a non-thermal fixed point of the evolution. The 
approach of this state can be identified by the build-up of a characteristic power-law single-particle 
spectrum n(k) ~ k~ 5 . Alternatively, for a weak flux in energy space, during the condensation process 
phase defects reminiscent of superfluid turbulence appear but a separation of the incompressible and 
compressible components, as well as the k~ 5 scaling of the low-energy modes are absent. 
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Introduction. In view of the rapid progress in exper- 
imental control of Bose-Einstein condensates (BEC) in 
dilute atomic gases, a growing community is investigat- 
ing the nonequilibrium dynamics of these systems. One 
of the most fundamental and long debated questions con- 
cerns the formation of a BEC from a completely disor- 
dered initial state [1-20]. A particularly interesting as- 
pect is the role of superfluid turbulence in this process [3- 
5, 14, 16]. Understanding the possible different paths to 
a BEC is of fundamental interest way beyond the realm 
of ultracold gases, from the phenomenology of the solid 
state up to the highest energies, e.g., in heavy-ion colli- 
sions [21-25] or early-universe evolution [21, 26-29]. 

Bose-Einstein condensation in a non-equilibrated and 
under-cooled gas can have the characteristics of a turbu- 
lent inverse cascade [3, 4, 14], corresponding to a quasi- 
local transport process in momentum space, into the low- 
energy modes of the Bose field. At energies below the 
regime where this cascade is expected to be described by 
kinetic weak- wave-turbulence theory [3, 4], a period with 
superfluid turbulence can appear which manifests itself 
through the appearance of vortex tangles [5, 14]. 

In this article, we demonstrate that the superfluid 
turbulence period can appear in two distinctly different 
forms. Suppose we start with a dilute, incoherent homo- 
geneous gas in three dimensions, at a temperature close 
to and above the BEC phase transition. Apply a cool- 
ing quench by removing particles from the higher-energy 
modes. After this, the system re-equilibrates to thermal 
equilibrium, with a temperature given by the total en- 
ergy of the system after the quench which we assume to 
be below the BEC critical temperature. 

Two different scenarios are possible: If a sufficiently 
small amount of energy is removed, the induced scatter- 
ing of particles into the low-energy modes builds up a 
thermal Rayleigh- Jeans distribution in a quasi-adiabatic 
way. The chemical potential increases, and a fraction 
of particles is deposited in the lowest mode, forming a 




FIG. 1: (Color online) Depending on the strength a of an ini- 
tial cooling quench, the gas can thermalize directly to a Bose- 
Einstein condensate, or it can first approach and critically 
slow down near a non-thermal fixed point (NTFP). There it 
is characterized by a self-similar particle spectrum n(k) ~ k~ 5 
and strong Vinen superfluid turbulence. 

BEC. During this process, tangles of defect lines can be 
found in the Bose field by filtering out short-wavelength 
fluctuations [14]. While the field has a phase winding 
around these defects, the flow pattern in their vicinity is 
distorted by strong phase fluctuations. 

In the second scenario, after a sufficiently strong cool- 
ing quench, the resulting super-thermal population of the 
higher-energy modes results in a vigorous transport to- 
wards lower energies that has the form of a strong- wave- 
turbulence inverse cascade. This cascade induces a long- 
lived, power-law single-particle spectrum n(k) ~ 
with an exponent £ = 5 distinctly larger than the expo- 
nent £ = 2 which characterizes a thermal Rayleigh- Jeans 
distribution. The emergence of the strong cascade can 
be explained by the dominance of incompressible, trans- 
verse, vortical superfluid flow over compressible, longitu- 
dinal sound excitations and density fluctuations in the 
respective regime of wave lengths. The power law k~ 5 
is traced back to the flow pattern around the vortex 
lines [33] in the Vinen tangles [30] which become visi- 
ble without filtering out short- wavelength fluctuations. 
The possibility of a strong cascade was pointed out in 
Refs. [27, 28, 31], where this quasi-stationary scaling solu- 
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FIG. 2: (Color online) Build-up of the condensate fraction 
n(k = 0,t)/N on a double- logarithmic scale, for different 
strengths of the initial cooling quench, parametrized by the 
exponent a, cf. Eq. (2). Depending on a, different power- law 
behavior is seen. Standard averaging errors are shown. 



tion was referred to as a non-thermal fixed point (NTFP). 
The two possible paths to BE condensation are shown 
schematically in Fig. 1. Whether the system, during the 
condensation dynamics, can approach the NTFP where 
the process critically slows down, or whether it moves in 
a direct way to thermal equilibrium depends on the initial 
conditions, i.e., on the strength of the cooling quench. We 
refer to the condensation process which takes the 'detour' 
via the NTFP as hydrodynamic BE condensation because 
of the dynamical separation of incompressible and com- 
pressible components in momentum space, in contrast to 
the direct process where this separation is absent. 

Semiclassical simulations. To reveal this dynamics we 
study a dilute Bose gas, in the classical- wave limit, using 
the Gross-Pitaevskii equation (GPE) (ft = 1) 



i<9 t 0(x,£) 
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We consider a gas of TV atoms in a box of size L 3 , 
with periodic boundary conditions and mean density 
n = N/L 3 . We measure length in units of the healing 
length £ = (2mgn)~ 1 ^ 2 and time in units of r = ra£ 2 . 
Simulations were done on a cubic grid with 256 3 points. 

The initial field in momentum space, 0(k, 0) = 
^ n(k, 0)exp{i(/?(k, 0)}, is parametrized in terms of a ran- 
domly chosen phase y?(k, 0) G [0, 2tt) and a density 
n(k, 0) = /(fc)z/k, with > drawn from an expo- 
nential distribution -P(^k) = exp(— z/jj f° r each k. We 
choose the structure function 



/(*) 



fa 



(2) 



for different values of a, with cutoff (k £) a = 0.2/0.44 a 
and normalization f a = 400/0. 44 a . We compare results 
for a range of different cooling quenches defined by the 
power-laws a = 2.5, . . . , 10.0, varying the total number 
between N = 10 9 (a = 2.5) and N = 4.3 x 10 8 (a = 10). 




o 



Radial momentum k£ 

FIG. 3: (Color online) Particle momentum spectra at four 
different times, on a double-logarithmic scale as functions of 
the radial momentum k = |k| for different a = 2.5 (orange), 
3.0 (blue), 4.0 (red), 6.0 (grey), 10 (black), from top to bottom 
at k£ — 1. Averages over 3.. 7 runs. A steep infrared power- 
law n(k) ~ k~ 5 appears for a > 3. 



Evolution of the zero mode. The initial power-law fall- 
off of f(k) chosen above is close to or steeper than that 
expected by a self-similar solution of the wave Boltzmann 
equation, a ~ 2.4 [7], corresponding to an inverse parti- 
cle cascade in weak wave turbulence theory [3, 4]. a = 2 
would correspond to a stable initial thermal distribution 
at finite chemical potential, and thus any a > 2 is re- 
quired to set off a re-equilibration to a lower tempera- 
ture. In Fig. 2 we show the ensuing time evolution of the 
condensate occupation number 77,(0, t) for the different a. 
In each case, the evolution leads to a BEC characterized 
by a non- vanishing ratio n(0,t)/N. As we keep f a and k$ 
constant n(0,i)/N grows with a. This is because larger 
a cut off more high-momentum particles and leave less 
energy to be thermally redistributed. Power-law growth 
n(0, t) rsj t [3], ~ t 3 [6], and ~ t 2 is seen. 

Momentum distribution. Beyond the zero mode, the 
occupation spectrum of the non-zero momentum modes 
allows to follow the dynamical process of potentially tur- 
bulent particle transport to lower energies. In Fig. 3 we 
show the time evolution of the single-particle distribution 
n(k,i) over the radial momenta k = |k| at four different 
times and for the different a. During the initial evolu- 
tion (t < 10 2 r) the mode occupations gradually spread to 
lower wave numbers, at the same time depositing energy 
into the high-momentum tail. Our dynamic range does 
not allow to identify a weak-turbulence inverse particle 
cascade with n(k) ~ k~ 2A [7]. At late times, the spec- 
tra developing from the different initial a differ strongly. 
For a > 3, the distribution develops a bimodal structure, 
with a power-law behavior n(k) ~ k~ 5 in the infrared 
(IR) and n(k) ~ k~ 2 in the UV. At very long times, this 
bimodal structure decays towards a global n(k) ~ k~ 2 
(not shown). For a < 3, the distribution directly reaches 
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FIG. 4: (Color online) Decomposition of particle spectra n(k) 
(black) into incompressible (red), compressible (blue), and 
'quantum-pressure' (grey) components, see text for defini- 
tions, on a double-log scale, at four different times, (a) Weak 
cooling quench, a = 2.5. Average of 3 runs (b) Strong quench, 
a = 10. Average of 7 runs 

a thermal Rayleigh- Jeans scaling n(k) ~ T/k 2 . We point 
out that the strongly nonthermal k~ 5 power law was 
recently predicted in a non- perturb at ive quantum-field- 
theory approach to strong wave turbulence [21, 31]. Its 
origin has been related to the presence of a dilute gas of 
vortices [32, 33]. 

Hydrodynamic condensation. To interpret our results 
in the context of superfluid turbulence we analyse kinetic- 
energy spectra as in [34]. The kinetic energy E^ n = 
/ dx(|V</>(x,f)| 2 )/2 is split, E kin = E v + E q , into the 
'classical' component E v = J dx (| y / nv| 2 )/2 with velocity 
v = Vy?/m, and a 'quantum-pressure' component E q = 
j dx (|Vy / n| 2 )/2. Radial particle spectra 

n 5 {k)= l - J dtt(|w,(k)| 2 ), 5 = v,q. (3) 

in terms of the generalised velocities w v = a/tiv and 
w q = V^Jn are furthermore decomposed by w v = Wi + 
w c into 'incompressible' (V • Wi = 0) and 'compressible' 
(V x w c = 0) parts to distinguish vortical superfluid and 
sound excitations of the gas, respectively. In Fig. 4a, 
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FIG. 5: (Color online) Evolution of the fraction A of inte- 
grated compressible to incompressible components below a 
momentum scale k\^ — 0.35, for different initial conditions. 
Note the semi-log scale. For a < 4, A(£) stays approximately 
constant. For a > 4, A(t) approaches zero at intermediate 
times, signaling the non-thermal fixed point and a superfluid 
hydrodynamic condensation process. Averages over 1...7 runs. 
Inset: A(t = 606r) as a function of a. 

we show the evolution of the different components n$(k), 
5 G {i, c, q}, for the weak initial quench a — 2.5. At early 
times, t < 10 3 r, due to the absence of phase coherence 
[33] the resulting spectra do not add up to the single- 
particle spectrum n(k) ^ n\{k) + n c (k) + n q (k). Com- 
pressible and incompressible components are roughly of 
equal magnitude while the quantum pressure component 
is insignificant on all scales. n(k) grows in the regime 
of low momenta while phase coherence is being estab- 
lished and a condensate fraction appears (Fig. 2). For 
the case of the strongly non-thermal initial distribution, 
a = 10, the evolution is shown in Fig. 4b. In contrast to 
Fig. 4a, two macroscopic flows can be observed, one to 
the UV, and one to the IR. Conservation of particle and 
energy imply immediately that, when sent out from the 
regime of intermediate frequencies, energy is deposited in 
the UV while particles are predominantly transferred to 
the IR. This leads to an inverse particle cascade with ap- 
proximately /^-independent radial particle flux Q(k) = Q 
and a corresponding direct energy cascade to the UV 
[33]. The inverse particle cascade reflects strong wave 
turbulence, characterized by n(k) ~ k~ 5 [31]. The de- 
composition in Fig. 4a makes clear that this power-law is 
caused by incompressible excitations only [32, 33], estab- 
lishing a dominantly ideal hydrodynamic (superfluid) BE 
condensation process. In the UV, the excitations follow 
a thermal n(k) ~ k~ 2 and are dominated by the com- 
pressible and quantum-pressure components. 

Our results show that during the hydrodynamic con- 
densation process incompressible flow temporally domi- 
nates at the expense of compressible excitations in the IR 
regime. The opposite occurs for the compressible exci- 
tations in the UV. This dynamical separation marks the 
approach to the NTFP and disappears when the system 
finally approaches thermal equilibrium. The evolution of 
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FIG. 6: (Color online) Vortex tangle structures emerging in 
the gas for the two different extremes of initial conditions, 
a = 2.5, at time t = 2424r (left) and a = 10, at time t = 606r 
(right). 



the integrated fraction A(t) — N c (k\,t)/Ni(k\,t), with 
Ns(k\,t) = i| k | <fcA dkn^(k, t), of compressible to incom- 
pressible occupations below a momentum scale k\^ = 
0.35 is shown in Fig. 5. Initially, A ~ 0.5. Starting from 
a < 3, A stays approximately constant while for a > 3, 
A decreases for t > 10 2 r towards zero before increasing 
again when thermal equilibrium is approached. The in- 
set shows A(t = 606) as a function of a in which one 
identifies a transition to a separation of the components, 
depending on the strength of the initial cooling quench. 

Defect formation. In Fig. 6 we show the three dimen- 
sional distribution of points where the density falls below 
0.2% of the average density n, for the systems quenched 
with a = 2.5 and a = 10, at the times t — 2A2At and 
606r, respectively. We filtered out modes with wavenum- 
ber larger than kt; = 0.45 but in the hydrodynamic case 
the high-momentum fluctuations barely distort the fig- 
ure. The vortex tangles corroborate the findings of [14] 
for both cases of a. However, a remarkable difference ex- 
ists in the distribution of the phase angle y?(x) of the Bose 
field as can be inferred from Fig. 4. While in the direct 
condensation case, the circular flow has strong longitu- 
dinal (compressible) fluctuations, in the hydrodynamic 
condensation via the NTFP macroscopic quantized vor- 
tical flow and thus Vinen turbulence of the superfluid are 
built up. Fig. 2 indicates a change in the power-law rise 
n(0, t) ~ t v from v ~ 3 to v ~ 2 to be associated with 
the approach of the NTFP. Similar features have been 
found to characterize the dynamics of a two-dimensional 
superfluid near a NTFP [35]. 

Summary. We find Bose-Einstein condensation in a 
dilute ultracold gas to occur in two different forms. Hy- 
drodynamic condensation involves a distinct dominance 
of incompressible flow bearing chaotic (Vinen) turbu- 
lence, corresponding to an approach of a non-thermal 
fixed point. In this regime, particles can not be deposited 
quickly enough into the zero mode and form an excess 
population with a characteristic power-law fall-off within 
the low-energy modes. In contrast, direct condensation 
can exhibit the appearance of vortical motion but lacks 



its separation from compressible sound-like excitations. 
The found characteristic power-law growth of the con- 
densate fraction allows to distinguish the regimes in ex- 
periment. 

Acknowledgements. We thank N. Berloff, J. Berges, 
J. P. Blaizot, S. Erne, M. Karl, L. McLerran, N. Philipp, 
J. Schole, D. Sexty, and B. Svistunov for useful discus- 
sions. Work supported by Deutsche Forschungsgemein- 
schaft (GA677/7,8), University of Heidelberg (FRON- 
TIER), and Helmholtz Association (HA216/EMMI). 



[1] E. Levich and V. Yakhot, Phys. Rev. B 15, 243 (1977); 

J. Phys. A: Math. Gen. 11, 2237 (1978). 
[2] H. T. C. Stoof, Phys. Rev. Lett. 66, 3148 (1991). 
[3] B. Svistunov, J. Mosc. Phys. Soc. 1, 373 (1991). 
[4] Y. Kagan, B. V. Svistunov, and G. V. Shlyapnikov, [Zh. 

Eksp. Teor. Fiz. 101, 528 (1992)] Sov. Phys. JETP 74, 

279 (1992). 

[5] Y. Kagan and B. V. Svistunov, [Zh. Eksp. Teor. Fiz. 105, 
353 (1994)] Sov. Phys. JETP 78, 187 (1994). 

[6] K. Damle, S. Majumdar, and S. Sachdev, Phys. Rev. A 
54, 5037 (1996). 

[7] D. Semikoz and I. Tkachev, Phys. Rev. D 55, 489 (1997). 

[8] Y. Kagan and B. V. Svistunov, Phys. Rev. Lett. 79, 3331 
(1997). 

[9] D. Snoke and J. Wolfe, Phys. Rev. B 39, 4030 (1989). 
[10] C. W. Gardiner, M. D. Lee, R. J. Ballagh, M. J. Davis, 

and P. Zoller, Phys. Rev. Lett. 81, 5266 (1998). 
[11] H. Miesner, et al., Science 279, 1005 (1998). 
[12] P. Drummond and J. Corney, Phys. Rev. A 60, 2661 

(1999). 

[13] M. Kohl, et al., Phys. Rev. Lett. 88, 080402 (2002). 
[14] N. G. Berloff and B. V. Svistunov, Phys. Rev. A 66, 
013603 (2002). 

[15] C. Connaughton, et al., Phys. Rev. Lett. 95, 263901 
(2005). 

[16] S. Nazarenko and M. Onorato, Physica D: Nonlin. Phen. 
219, 1 (2006). 

[17] S. Ritter, et al., Phys. Rev. Lett. 98, 090402 (2007). 

[18] C. N. Weiler, et al., Nature 455, 948 (2008). 

[19] E. V. Kozik and B. V. Svistunov, J. Low Temp. Phys. 
156, 215 (2009). 

[20] R. P. Smith, et al., arXiv: 1112.4457 [cond-mat. quant- 
gas] (2011). 

[21] J. Berges, S. Schemer, and D. Sexty, Phys. Lett. B681, 
362 (2009). 

[22] M. Carrington and A. Rebhan, Eur. Phys. J. C71, 1787 

(2011) . 

[23] K. Fukushima and F. Gelis, Nucl. Phys. A874, 108 

(2012) . 

[24] J.-P. Blaizot, et al., Nucl. Phys. A873, 68 (2012). 
[25] J. Berges and D. Sexty, Phys. Rev. Lett. 108, 161601 
(2012). 

[26] R. Micha and I. I. Tkachev, Phys. Rev. Lett. 90, 121301 
(2003). 

[27] J. Berges, A. Rothkopf, and J. Schmidt, Phys. Rev. Lett. 

101, 041603 (2008). 
[28] J. Berges and G. Hoffmeister, Nucl. Phys. B813, 383 

(2009). 



5 



[29] T. Gasenzer, B. Nowak, and D. Sexty, Phys. Lett. B710, 
500 (2012). 

[30] G. E. Volovik, J. Low Temp. Phys. 136, 309 (2004). 
[31] C. Scheppach, J. Berges, and T. Gasenzer, Phys. Rev. A 

81, 033611 (2010). 
[32] B. Nowak, D. Sexty, and T. Gasenzer, Phys. Rev. B 84, 

020506(R) (2011). 



[33] B. Nowak, J. Schole, D. Sexty, and T. Gasenzer, Phys. 

Rev. A 85, 043627 (2012). 
[34] C. Nore, M. Abid, and M. E. Brachet, Phys. Rev. Lett. 

78, 3896 (1997). 
[35] J. Schole, B. Nowak, and T. Gasenzer, arXiv: 1204.2487 

[cond-mat. quant-gas] (2012). 



